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We show that large amplitude, coherent acoustic phonon wavepackets can be generated and 
detected in In^Gai-ajN/GaN epilayers and heterostructures in femtosecond pump-probe differential 
reflectivity experiments. The amplitude of the coherent phonon increases with increasing Indium 
fraction x and unlike other coherent phonon oscillations, both amplitude and period are strong 
functions of the laser probe energy. The amplitude of the oscillation is substantially and almost 
instantaneously reduced when the wavepacket reaches a GaN-sapphire interface below the surface 
indicating that the phonon wavepackets are useful for imaging below the surface. A theoretical 
model is proposed which fits the experiments well and helps to deduce the strength of the phonon 
wavepackets. Our model shows that localized coherent phonon wavepackets are generated by the 
femtosecond pump laser in the epilayer near the surface. The wavepackets then propagate through a 
GaN layer changing the local index of refraction, primarily through the Franz-Keldysh effect, and as 
a result, modulate the reflectivity of the probe beam. Our model correctly predicts the experimental 
dependence on probe-wavelength as well as epilayer thickness. 



I. INTRODUCTION 

Heterostructures of GaN and InGaN are important 
materials owing to their applications to blue laser diodes 
and high power electronics^ Strong coherent acous- 
tic phonon oscillations have recently been detected in 
InGaN/GaN multiple quantum wells* 2 ^ These phonon 
oscillations were much stronger than folded acoustic 
phonon oscillations observed in other semiconductor 
superlatticesA^ 

InGaN/GaN heterostructures are highly strained at 
high In concentrations giving rise to large built-in piezo- 
electric fields and the large strength of the coher- 
ent acoustic phonon oscillations was attributed to the 
large strain and piezoelectric fields £ 

In this paper, we report the generation of strong lo- 
calized coherent phonon wavepackets in strained layer 
In x Gai_ x N/GaN epilayers and heterostructures grown 
on GaN and Sapphire substratesiii By focusing high rep- 
etition rate, frequency-doubled femtosecond Ti:Sapphire 
laser pulses onto strongly strained InGaN/GaN het- 
erostructures, we can, through the carrier-phonon in- 
teraction, generate coherent phonon wavepackets which 
are initially localized near the epilayer/surface but then 
propagate away from the surface/epilayer and through a 
GaN layer. As the wavepackets propagate, they modu- 
late the local index of refraction and can be observed in 
the time-dependent differential reflectivity of the probe 
pulse. There is a sudden drop in the amplitude of the re- 
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flectivity oscillation of the probe pulse when the phonon 
wavepacket reaches a GaN-sapphire interface below the 
surface. Theoretical calculations as well as experimental 
evidences support this picture; the sudden drop of ampli- 
tude when the wave encounters the GaN-sapphire inter- 
face cannot be explained if the wavepacket had large spa- 
tial extent. When the wavepacket encounters the GaN- 
sapphire interface, part of the wave gets reflected while 
most of it gets transmitted into the sapphire substrate, 
depending on the interface properties and the excess en- 
ergy of the exciting photons. This experiment illustrates 
a non-destructive way of generating high pressure ten- 
sile waves in strained heterostructures and using them to 
probe semiconductor structure below the surface of the 
sample. Since the strength of this non-destructive wave 
is determined by the strain between GaN and InGaN, 
it is likely that even stronger coherent phonons can be 
generated in InGaN/GaN digital alloys grown on a GaN 
substrate. 



II. EXPERIMENTAL RESULTS 

In the experiments, frequency-doubled pulses of mode- 
locked Tksapphire lasers are used to perform reflec- 
tive pump-probe measurements on four different sam- 
ple types which are shown in Fig. Type (I) InGaN 
epilayers; Type (II) InGaN/GaN double quantum wells 
(DQWs); Type (III) InGaN/GaN single quantum well 
(SQW); and Type (IV) InGaN/GaN light-emitting diode 
(LED) structures. The peak pump power is estimated 
to be 400 MW/cm 2 , corresponding to a carrier density 
of 10 19 cm~ 3 and the doubled pulse width is 250 fs. 
All samples were grown on a c-plane sapphire substrate 
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FIG. 1: The simple diagram of the sample structures used in these experiments. (a)InGaN epilayer (type I), (b)InGaN/GaN 
double quantum wells (type II), (c)GaN single quantum well (type III), (d)InGaN/GaN blue light-emitting diode structure 
(type IV). 



by metal organic chemical vapor deposition. The In- 
GaN epilayers shown in Fig. da) consist of 1 /im GaN 
grown on a sapphire substrate and capped with 30 nm 
of In^Gai-^N with In composition, x, varying from 0.04 
to 0.12. The DQW sample shown in Fig. ^b) consists 
of GaN (1 /im), double quantum wells of Ino.12Gao.ssN 
(1-16 nm), a barrier of GaN (1-16 nm), and a GaN cap 
layer (0.1 /im). The SQW sample in Fig. QJc) consists 
of GaN (2 /im), an Ino.12Gao.ssN well (24 nm), and a 
GaN cap layer (0.1 /im). The blue LED structure shown 
in Fig.^d) consists of n-GaN (4 /im), 5 quantum wells 
of Ino.15Gao.85N (2 nm) and 4 barriers of GaN (10 nm), 
p-Alo.1Gao.9N (20 nm), andp-GaN (0.2 /im). 

Differential reflection pump-probe measurements for 
the In^Gai-^N epilayers are shown in Fig. [TJa). Fig- 
ure [21 shows the oscillatory component of the measured 
probe differential reflectivity for the InGaN epilayers 
(type I) with various In concentrations, x. For com- 
parison purposes, we performed differential reflectivity 
measurements on a pure GaN HVPE grown sample in 
order to show that no differential reflectivity oscillations 
are present in the absence of strain and an epilayer. The 
energy of the pump laser was varied between 3.22 and 
3.35 eV, to keep the excess carrier energy above the In- 
GaN band gap but below the GaN band gap. We note 
that if the laser energy was below the InGaN band gap, 
no signal was detected. Therefore, carrier generation is 
essential to observing the oscillations, unlike a recent co- 
herent optical phonon experiment in GaN—. The inset 
shows the pump-probe signal prior to the background 
subtraction for x = 0.12. The background results from 
the relaxation of the photoexcited electrons and holes. 
The oscillations are quite large, on the order of 10 -2 - 
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FIG. 2: The oscillatory component of the differential reflec- 
tion pump-probe data for the In^Gai-^N epilayers with var- 
ious In composition (x=0.04, 0.08, 0.10, and 0.12). For com- 
parison, differential reflection in a pure GaN HPVE grown 
sample is shown. The reflection signal prior to the background 
subtraction for x — 0.12 is shown in the inset. 
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FIG. 3: The oscillation traces of a SQW (III) at different 
probe energies. The pump energy is centered at 3.26 eV. The 
bottom curve has been magnified 30 times. The inset shows 
the oscillation amplitude as a function of probe photon energy 
on a logarithmic scale. 



10~ 3 and the period is 8-9 ps, independent of the In 
composition but dependent on the probe photon energy. 
The amplitude of the oscillation is approximately propor- 
tional to the In concentration indicating that the strain at 
the InGaN/GaN interface is important. The observed pe- 
riod is approximately r = X/2C s n, where A is the probe 
beam wavelength, C s the longitudinal acoustic sound ve- 
locity, and n the refractive index of GaN 13 . 

We performed two-color pump-probe experiments for a 
type III InGaN SQW sample as shown in Fig.QIb). Fig.f5| 
shows the differential reflectivity oscillations for different 
probe energies. Note that the period of the oscillation 
changes and is proportional to the probe wavelength. In 
addition, the amplitude of the differential reflectivity os- 
cillation decreases as the detuning (with respect to the 
pump) becomes larger. The inset shows the oscillation 
amplitude as a function of the probe energy in a logarith- 
mic scale and there is an ~2-order-of magnitude decrease 
in differential reflectivity when the probe energy changes 
from 3.26 eV to 1.63 eV. 

Interestingly, Fig. 0] shows that the amplitude of the 
oscillatory component of the differential reflectivity is in- 
dependent of the bias voltage, even though the carrier 
lifetime changes dramatically with voltage bias. Fig. fj] 
shows the bias dependent acoustic phonon differential re- 
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FIG. 4: (a)Pump-probe differential reflectivity for the blue 
LED structure (type IV) varying the external bias at a pump 
energy of 3.17 eV. The decay time of the background signal is 
drastically reduced as the bias increases, (b) The oscillatory 
amplitude does not change much with bias. 



flectivity oscillations in a type IV blue LED structure (see 
Fig-IHd)) at a pump energy of 3.17 eV. The lifetime of 
the background signal drastically decreases as the bias 
increases as shown in Fig.@Ja). This is due to the carrier 
recombination time and the decrease in the tunneling es- 
cape time in the strong external bias regime^*. On the 
other hand, the amplitude and frequency of the oscil- 
latory component of the differential reflectivity doesn't 
change much with bias voltage [Fig.|3Jb)]. Since the ob- 
served reflectivity oscillation is independent of the carrier 
lifetime for lifetimes as short as 1 ps due to ultrafast tun- 
neling [bottom curve of Fig. Ufa)] , it implies that once 
the source that modulates the experimentally observed 
reflectivity is launched by the sub-picosecond generation 
of carriers, the remaining carriers do little to affect the 
source. This suggests that the reflectivity oscillation is 
due to the strain pulse which is generated at short times 
once the pump excites the carriers and modulates the 
lattice constant. 

Figure [S] shows the well- and barrier-dependent acous- 
tic phonon differential reflectivity oscillations in the type 
II DQW samples (see Fig.^b)) at a pump energy of 3.22 
eV. The amplitude increases as the well width increases. 
However, the oscillation amplitude of the differential re- 
flectivity doesn't change much with the barrier width. 
This means that the generation of the acoustic phonons 
is due to the InGaN well and not the GaN barrier. This 
also verifies that the oscillation is due to the strain in 
InGaN layer. 

Interesting results are seen in the long time behavior 
of the reflectivity oscillations shown in Fig. OJa). The 
long time scale reflectivity oscillation is plotted for the 
epilayer (I), DQW (II), SQW (III), and the blue LED 
structure (IV) at 3.29 eV (below the GaN band gap). As- 
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FIG. 5: The oscillatory component of the pump-probe differ- 
ential transmission traces of DQW's (II) at 3.22 eV. The left 
figure shows the well width dependence and the right figure 
shows the barrier width dependence. 



tonishingly, the oscillation amplitude abruptly decreases 
within one cycle of an oscillation at a critical time which 
appears to scale with the thickness of the GaN layer in 
each sample. In addition, the slope of the GaN thick- 
ness vs. the critical time is very close to the known 
value of the sound velocity in GaN [inset of Fig. EOA)} 2 . 
Fig. E{b) shows results when the probe laser energy is 
changed to 3.44 eV which is above the GaN band gap. 
Then the laser probe is sensitive to coherent phonon os- 
cillations only within an absorption depth of the surface. 
We see that the amplitude of the reflectivity oscillation 
exponentially decays with a decay time of 24.2 ps corre- 
sponding to a penetration depth in GaN of about 0.17 
micron (=24.2 psx7000 m/s). The oscillation reappears 
at 260 ps for epilayer (I) and 340 ps for DQW (II). This 
is twice the critical time for the oscillations to disap- 
pear when the photon energy is 3.29 eV. This further 
shows that the probe pulse is sensitive to the coherent 
acoustic wave. The "echo" in the probe signal results 
from the partial reflection of the coherent phonon off the 
GaN/sapphire interface. 
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FIG. 6: (a)The long time-scale oscillation traces of an epi- 
layer (I), a DQW (II), a SQW (III), and blue LED structure 
(IV) . The inset shows the GaN thickness between the sapphire 
substrate and InGaN active layer of the sample as a function 
of the die out time of the oscillations. The solid line indicates 
that the velocity of the wavepacket in the GaN medium is 
about 7000 m/s. (b) The oscillation traces of epilayer and 
DQW [top two curves in (a)]at 3.44 eV which corresponds to 
a probe laser energy above the band gap of GaN 



III. THEORY 

To explain the experimental results discussed in the 
last section, we have developed a theoretical model of 
the generation, propagation and detection of coherent 
acoustic phonons in strained GaN/InGaN heterostruc- 
tures. The pump laser pulse generates a strain field that 
propagates through the sample which, in turn, causes a 
spatio-temporal change in the index of refraction. This 
change is responsible for the oscillatory behavior seen 
in the probe-field reflectivity in various semiconductor 
heterostructures. An approximate method of solving 
Maxwell's equations in the presence of spatio-temporal 
disturbances in the optical properties and obtaining the 
reflectivity of the probe field in thin films excited by pi- 
cosecond pump pulses can be found in Thompsen et. al^ 

The spatio-temporal disturbance of the refractive in- 
dex is caused by the propagating coherent phonon 
wavepackets. Thus, an essential ingredient in under- 
standing the probe reflectivity is a model for the gen- 
eration and propagation of the very short strain pulse 
in the sample. Recently, a microscopic theory explain- 
ing the generation and propagation of such strain pulses 
was reported by Sanders et. alpi£*i& where it was shown 
that propagating coherent acoustic phonon wavepackets 
are created by the nonequilibrium carriers excited by the 
ultrafast pump pulse. The acoustic phonon oscillations 
arise through the electron-phonon interaction with the 
photoexcited carriers. Both acoustic deformation poten- 
tial and piezoelectric scattering were considered in the 
microscopic model. It was found that under typical ex- 
perimental conditions, the microscopic theory could be 
simplified and mapped onto a loaded string model. Here, 
we use the string model of coherent phonon pulse gener- 
ation to obtain the strain field seen by the probe pulse. 
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First, we solve Maxwell's equations to obtain the 
probe reflectivity in the presence of a generalized spatio- 
temporal disturbance of the index of refraction. Let rib 
be the index of refraction without the strain which is real 
because initially the absorption can be neglected and let 
5n = 5n + in, be the propagating change in the index 
of refraction due to the strain. When the effect of the 
change of the index of refraction is taken into account, 
the probe field with energy u> can be described by the 
following generalized wave equation 

d 2 E(z,t) oj 



dz 2 



+ ^-[n h + 5n(z,t)] 2 E(z,t) = (1) 



where E(z, t) is the probe field in the slowly varying 
envelope function approximation and lj is the central 
frequency of the probe pulse. Eqn. Q is obtained 
from Maxwell's equations assuming that the polariza- 
tion response is instantaneous and that the probe pulse 
obeys the slowly varying envelope function approxima- 
tion. Since \5n\ <C «b under typical conditions, Eqn. 
can be cast into 



d 2 E(z,t) 
dz 2 



where k = uj/c is the probe wavevector. To relate the 
change of the index of refraction to the strain field, 
rj(z,t) = dU(z,t)/dz, we assume \Sn\ <C rih and adopt 
the linear approximation^ 



. . dn 
on{z,t) = —<n(z,t). 



(3) 



We view Eqn. J2} as an inhomogencous Hclmholtz 
equation and obtain the solution using the Green's func- 
tion technique. The desired Green's function is deter- 
mined by solving 



d 2 G(z,z') 
dz 2 

and the result is 



+ (n h k) 2 G(z,z r ) = 8{z-z') 



(4) 



G(z,z') = -- — -exp(in h k\z - z'\) . (5) 
ln\ } k 



+ (nbfc) E(z,t) = — 2ribfc 5n(z,t)E(z,t) (2) Then, the solution to Eqn. 10) can be written as 



E(z,t) = E h {z,t)+ / dz'G(z,z') {~2n b k 2 Sn{z',t) E(z',t)} 

J — OO 

dz'G{z,z'){~2n h k 2 6h(z',t) E h (z',t)} 

-OO 

I 



(6) 
(7) 



where we have chosen the lowest Born series in the last 
line. In Eqn. J7J), E^ is the homogeneous solution which 
takes the form Eh(z,t) — E^(z,t) exp {i(n\^kz — Lot)} 
where E^(z, t) is a slowly varying envelope function. This 
describes the probe pulse moving to the right in the sam- 
ple without optical distortion. 

We now apply this approximate solution to our struc- 
ture where the interface between air and the sample is 
chosen at z = 0. In the air, where z < 0, there is an 
incident probe pulse traveling toward the sample as well 



as a reflected pulse. The electric field in the air, E < (z, t), 
can thus be written as the sum 

E<{z,t) = Ei(z,t) e l( - kz -^ +E T (z,t) e ~^ kz+ ^ (8) 

where Ei(z,t) and E x {z,t) are the slowly varying enve- 
lope functions of the incident and reflected probe fields, 
respectively. Inside the sample, z > 0, the solution is 
given as 



J 



E > (z,t) = E h (z,t)+ / dz' exp (iribk(z' — z)) ik 5h(z' , t) E h (z' , t) 
Jo 

= E t {z 7 t)e l( - nhkz -^ + jjT dz'exp(2in h kz')ik5fi(z\t)^E t (z,t)e' l{nbkz+ut) (9) 
= E t { Zl t)e l{nhkz - ut) +A{n h k 1 t)E t {z,t)e- l{nbkz+ut) (10) 



where we used Eh(z, t) = E t (z, t) exp(i(nbA:z — tot)) in the second step, assuming that E t (z, t) is nearly constant 
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within the slowly varying envelope function approxima- 
tion, and we define a reflected amplitude function 

/>oo 

A(nt,k,t) = / dz' exp [2in\Jcz') ikSn(z',t). (11) 
Jo 

The expression for the reflected amplitude function, 
A{n\,k, t), in Eqn. J^J says there is a frequency-dependent 
modulation of the amplitude in the reflected wave in the 
sample due to the propagating strain. 

Having determined the waves on both sides of the in- 
terface, we can now calculate the reflectivity. We apply 
the usual boundary conditions to the slowly varying en- 
velope functions and the results are written compactly as 



-1 l+A 
1 n h (l-A) 



E r 



Ei 



We solve this equation to obtain 



El 



rp+A 
l + r A 



r + A 



(12) 



(13) 



where tq = (1 — n,b)/(l + rib)- To the same order, we find 
that E t /Ei w to(l — TqA) where to — 2/ (l + n b ). It is now 
straightforward to calculate the differential reflectivity as 



AR _ \r + Ar\ 2 



Fol 



Fo 



—Re A (14) 



Finally, by substituting the linear law Eqn. J3J) into 
Eqn. lfTT|) and using Eqn. ((Til we get 

— =j o dzFfru)—— (15) 

where the sensitivity function, J-{z,uj), is defined as 



2k 



dn . . , dn . . 
— simzn^kz) + — cos(2nhkz) 
orj o-q 



(16) 



In Eqn. Ijl5(l , the differential reflectivity is expressed in 
terms of the lattice displacement, U(z,t), due to propa- 
gating coherent phonons. Sanders et. ali^l developed 
a microscopic theory showing that the coherent phonon 
lattice displacement satisfies a driven string equation, 



d 2 U(z,t) 
dt 2 



d 2 U(z,t) 
dz 2 



S(z,t), 



(17) 



where C s is the LA sound speed in the medium and 
S(z,t) is a driving term which depends on the photo- 
generated carrier density. The LA sound speed is related 
to the mass dens ity, p, and the elastic stiffness constant, 
C33, by C s = y/ C33 / p. The LA sound speed is taken 
to be different in the GaN/InGaN heterostructure and 
Sapphire substrate. For simplicity, we neglect the sound 
speed mismatch between the GaN and In x Gai_ x N layers. 



The driving function, 
given by 



S(z,t), is nonuniform and is 



(18) 



where the summation index, is, runs over carrier species, 
i.e., conduction electrons, heavy holes, light holes, and 
crystal field split holes, that are created by the pump 
pulse. Each carrier species makes a separate contribution 
to the driving function. The partial driving functions, 
S v (z,t), in piezoelectric wurtzite crystals depend on the 
density of the photoexcited carriers. Thus, 



a f j.\ , 1 J 9 47r|e| e 33 

S v (z,t) = ±- < a v — -\ 

P { oz 



p„{z,t), (19) 



where the plus sign is used for conduction electrons and 
the minus sign is used for holes. Here p u (z,t) is the pho- 
togenerated electron or hole number density, which is real 
and positive, p is the mass density, a„ are the deforma- 
tion potentials for each carrier species, 633 is the piezo- 
electric constant, and £oo is the high frequency dielectric 
constant. 

To be more specific, we will consider a SQW sample of 
the type shown in Fig. ^c) . We adopt a simple model 
for photogeneration of electrons and holes in which the 
photogenerated electron and hole number densities are 
proportional to the squared ground state particle in a box 
wavef unctions. The exact shape of the electron and hole 
number density profile is not critical in the present cal- 
culation since all that really matters is that the electrons 
and holes be strongly localized. The carriers are assumed 
to be instantaneously generated by the pump pulse and 
are localized in the In^Gai-^N quantum well. Having 
obtained a model expression for p v {z, t), it is straightfor- 
ward to obtain S(z,t) using Eqns. I|18fl and (|19fl . 

To obtain U(z,t), we solve driven string equations in 
the GaN epilayer and the Sapphire substrate, namely 



d 2 U(z,t) 



dt 2 



r 2 



d 2 U(z,t) 

dz 2 



and 



d 2 U{z,t) 
dt 2 



2 d 2 U(z,t) 
Ll dz 2 



S(z,t) (0<z<L) 
(20a) 

= (L < z < Z s ) 

(20b) 

where Co and C\ are LA sound speeds in the GaN and 
Sapphire substrate, respectively. In Eqn. (|20b|l . the Sap- 
phire substrate has finite thickness. To simulate coherent 
phonon propagation in an infinite Sapphire substrate, Z s 
in Eqn. I|20b|l is chosen large enough so that the propa- 
gating sound pulse generated in the GaN epilayer does 
not have sufficient time to reach z = Z s during the sim- 
ulation. If T sim is the duration of the simulation, this 
implies Z s > L + C\ T sim . 

Equations l|20a|) and (|20b|) are solved subject to initial 
and boundary conditions. The initial conditions are 



U(z,0) 



dU(z,0) 
dt 



0. 



(21) 
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At the GaN-air interface at z = 0, we assume the free 
surface boundary condition 



dU{0,t) 
dz 







(22a) 



since the air exerts no force on the GaN epilayer. The 
phonon displacement and the force per unit area are con- 
tinuous at the GaN-Sapphire interface so that 



U(L-e,t) = U{L + e,t) 



and 



Po C[ 



2 dU{L-e,t) _ ^ 2 8U(L + e,t) 

o ^ Pi c i ^ • 

oz dz 



(22b) 



(22c) 



The boundary condition at z — Z s can be chosen arbi- 
trarily since the propagating sound pulse never reaches 
this interface. For mathematical convenience, we choose 
the rigid boundary condition 



U(Z s ,t)=0. 



(22d) 



To obtain U(z,t) for general S(z,t), we first need to 
find the harmonic solutions in the absence of strain, i.e. 
S(z,t) = 0. The harmonic solutions are taken to be 

U n (z,t) = Wn(z) e-*"»* (w„>0) (23) 

and it is easy to show that the mode functions, W n (z), 
satisfy 



d 2 W n (z) t loI 

-o 



(h , +(%W n (z) = (0<z<L) (24a) 



and 



+ ^W n (z) = ( L < z < Z s ) (24b) 
az z C( 

Applying the boundary conditions 12211 we obtain the 
mode functions 



W n {z) 



with 



cos (u> n z/C ) 



if < z < L 



B n sin (uj n (Z s - z)/Ci) if L < z < Z s 

(25a) 



COS (uJnL/Cg) 

sin (uj n (Z s - L)/Ci) 



(25b) 



The mode frequencies, u> n , are solutions of the transcen- 
dental equation 



ro '' or 



PqCq 



PiC 



1 ^( "nfr-LY ^ (26) 



Ci 



which we solve numerically to obtain the mode frequen- 
cies, uj n (n = 0,1,2,...). The index, n, is equal to the 
number of nodes in the mode functions, W n (z). 



A general displacement can be expanded in terms of 
the harmonic modes as 



U(z,t) = J^ q n (t) W n (z 



(27) 



71=0 



Substituting equation (|27() for U(z,t) into (|20|l and tak- 
ing the initial conditions (|21|) into account, we find that 
the expansion coefficients, q n (t), satisfy a driven har- 
monic oscillator equation 



?„(*)+ u£ ?„(*)= On (*), 



(28) 



subject to the initial conditions q n (0) — <z«(0) = 0. The 
harmonic oscillator driving term Q n (t) is given by 



Qn(t) 



J Z 'dz W n (z) S(z,t) 
J Z 'dzW n (zr 



(29) 



In our simple displacive model for photogeneration of 
carriers, S(z,t) = S(z) 0(f) where O(i) is the Heaviside 
step function. In this case, the lattice displacement is 
explicitly given by 



U(z, t) = V -=■( 1 - cos(cj n t) ) W n (z) (30) 
with S n defined as 



Jo' dz W n (z) S(z) 
f Z 'dzW n (z)* 



(31) 



Using the lattice displacement H30(l . we obtain the time- 
dependent differential reflectivity at the probe frequency, 
u>, from equation (|15|) . The result is 



Ai? S 

-(w,f) = V -£( 1 - cosK t) ) R n {u) (32) 



R 



where 



Zs dz Hz , w) MM. (33) 
dz 



can be evaluated analytically. 

With the above formalism, we solve for the lattice 
displacement, U(z,t), for a coherent LA phonon pulse 
propagating in a multilayer structure consisting of a 
1.124 /im thick GaN epilayer grown on top of an in- 
finitely thick Sapphire substrate with the growth direc- 
tion along z. We take the origin to be at the GaN- 
air interface and the GaN-Sapphire interface is taken 
to be at z = L = 1.124 /xm. We assume that carriers 
are photogenerated in a single 240 A thick In^Gai-^N 
quantum well embedded in the GaN layer 0.1 /an be- 
low the GaN-air interface and 1 /im above the Sapphire 
substrate. Our structure thus resembles the SQW sam- 
ple shown in Fig. ^c). In the GaN epilayer, we take 
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tained from equation (|15|l . From Fig. da), the strains, 
dU(z,t)/dz, associated with the propagating pulses are 
highly localized and travel at the LA sound speed. Each 
pulse contributes a term to the differential reflectivity 
that goes like 



AR. . . u . 

—- (u>, t) ~ ^(Got, u) oc - sin 
R c 



Cot + cj)) (34) 



The period of the oscillations of T depends on the probe 
wavelength, A = 2itc/uj, with the result that the ob- 
served differential reflectivity oscillates in time with pe- 
riod, T — nc/(nbCnOj) = A/(2nbCo), where = 2.4 
is the index of refraction, and Co = 7857 m/s is the 
LA sound speed in GaN. For A = 377 nm (hu) = 3.29 
eV), this gives us T = 10 ps. The sensitivity function, 
J-(z,lu), defined in equation l|16|) is an oscillating func- 
tion in the GaN/InGaN epilayer and is assumed to vanish 
in the Sapphire substrate. Our computed differential re- 
flectivity is shown in Fig. |7Jb) for a probe wavelength of 
A = 377 nm. We find that the reflectivity abruptly atten- 
uates when the strain pulse enters the Sapphire substrate 
at t = 170 ps. The reflected strain pulses give rise to the 
weaker oscillations seen for t > 170 ps. These oscilla- 
tions are predicted to continue until the reflected pulses 
are again partially reflected off the Sapphire substrate at 
t = 430 ps. 



FIG. 7: Generation and propagation of coherent acoustic 
phonons photogenerated in a single InGaN well embedded 
in a free standing 1.124/im GaN epilayer grown on top of a 
Sapphire substrate. In (a) a contour plot of the strain field, 
dU(z,t)/dz, is shown as a function of depth below the GaN- 
air interface and the probe delay. In (b), the resulting differ- 
ential reflectivity induced by the strain field in (a) is shown 
as a function of the probe delay. 



C33 = 379 GPa and po = 6.139 gm/cm 3 — from which 
we obtain Co = 7857 m/s. For the Sapphire substrate, 
we take C33 = 500 GPa and pi = 3.986 gm/cm 3 i£ from 
which we find C\ = 11200 m/s. 

The results of our simulation are shown in Fig. \7\ 
A contour map of the strain, dU(z,t)/dz, is shown in 
Fig. [7Ja). We plot the strain as a function of the depth 
below the surface and the probe delay time. Photoexci- 
tation of electrons and holes in the InGaN quantum well 
generates two coherent LA sound pulses traveling in op- 
posite directions. The pulses are totally reflected off the 
GaN/air interface at z = and are partially reflected at 
the Sapphire substrate at z = 1.124 pm. Approximately 
95% of the pulse energy is transmitted and only 5% is 
reflected at the substrate. The speed of the LA phonon 
pulses is just the slope of the propagating wave trains 
seen in Fig. Efa.) and one can clearly see that the LA 
sound speed is greater in the Sapphire substrate. 

From the strain the differential reflectivity can be ob- 



IV. SIMPLE MODEL 

Since the coherent oscillation observed in the differ- 
ential reflectivity stems essentially from the strain pulse 
propagating into the layers, most phenomena can be un- 
derstood by a simple macroscopic model that is presented 
in this section. 

Instead of solving the loaded string equations for the 
strain to obtain a propagating disturbance in the refrac- 
tive index, the propagating strain pulse at a given mo- 
ment can be viewed as a thin strained layer in the sample, 
where the index of refraction is assumed to be slightly 
different from the rest of the sample. This situation is 
schematically depicted in Fig. [S] where a fictitious, thin 
GaN strained layer is located at z in the thick host GaN 
layer. The thickness of the strained layer, d, is approxi- 
mately the width of the traveling coherent phonon strain 
field, dU(z,t)/dz, and is to be determined from the mi- 
croscopic theory. From the last section, it was seen that 
the propagating strain field is strongly localized so that 
d is small. In the example of the last section, d is ap- 
proximately equal to the quantum well width. Here we 
assume the strain pulse has been already created near the 
air/GaN interface and do not consider its generation pro- 
cedure. We treat d as a phenomenological constant and 
also assume that the change in the index of refraction is 
constant. This strained GaN layer travels into the struc- 
ture with the speed of the acoustic phonon wavepacket 
Co = 7 x 10 3 m/s, so the location of the stained layer is 
given as z = Cqt where r is the pump-probe delay time. 
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FIG. 8: Propagating strained GaN layer in our simple model. 
The pump laser pulse creates a coherent acoustic phonon 
wavepacket in the InGaN layer near the air/GaN surface, 
which is modelled as a thin strained layer. The strained 
GaN layer propagates into the host GaN layer. The index 
of refraction in the strained layer is perturbed relative to the 
background GaN due to the strong strain induced piezoelec- 
tric field (Franz-Keldysh effect). 

Within the slowly varying envelope function approx- 
imation, the solutions to the Maxwell equation can be 
written as plane waves in each region, 

Ei(z, t) = a i e tt *"*- fa '* + bie- ikiZi - iut (35) 

where E{ is the electric field in the left layer of ith in- 
terface, and di and hi are the slowly varying amplitudes. 
The magnetic field is given by B. L ~ dEi/dz. For a nor- 
mally incident probe field, we apply the usual matching 
conditions on E and B to obtain 

(X )-*•(&)■ (36 » 

where the transfer matrix, Mi, is given explicitly as 
Mi = (37) 

1 / (1 | fc * +1 ) c »(fc'+i -ki)zi (1 *<+i ) c -»ffcM-i+k1gt \ 

2 I (1 - *<±i) e i(fc<+i+fc<)** (1 _|_ *j£tl.) e -»(fc<+i J ■ 

To apply this formula to our configuration, we normalize 
the incident amplitude to 1, let r be the reflected am- 
plitude in the air, and impose the boundary condition 
that there is only a transmitted wave into the Sapphire 
substrate with amplitude t and no reflected wave from 
the GaN-Sapphire interface back into the GaN epilayer. 
The latter assumption is reasonable since the microscopic 
theory of the previous section shows that only 5% of the 
pulse energy is reflected from the interface between the 
GaN and the Sapphire substrate. The total reflection 
and transmission amplitudes r and t for the GaN epi- 
layer structure are determined by 



0.02 
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FIG. 9: Differential reflection for different frequencies of the 
probe pulse. Plot (a) is experimental data and (b) comes from 
our calculations described in the text. 

where M = M1M2 • • • M n . The reflection amplitude is 
readily found to be 

r = M 2 \/Mu =r + Ar, (39) 

where ro is the background contribution without the 
strained layer. The total transmission amplitude is given 
as t = 1/Mn- We can now numerically calculate the 
differential reflectivity by substituting r into Eqn. 114|) . 

In Fig. we compare our theoretical model with exper- 
imental results. In Fig. IHJa), the oscillatory part of the 
differential reflectivity is plotted as a function of probe 
delay for three values of the probe photon energy. This 
is the same data that was shown in Fig. [3] where it was 
seen that the oscillation period of the coherent phonon re- 
flectivity oscillations are proportional to the probe wave- 
length. In Fig. IHIb), we plot the corresponding theo- 
retical differential reflectivities obtained from Eqn. H14[) . 
Comparing Figs.|§fa) andl^b), we see that our geomet- 
rical optics model successfully explains the observed rela- 
tion between the coherent phonon oscillation period and 
the probe wavelength. However, the calculated ampli- 
tudes of oscillation is inconsistent with the experimental 
data because in this case we have not taken into account 
of the change of index of refraction with respect to the 
probe wavelength which we will do later in the paper. 

In Fig.^5]we plot the differential reflectivity as a func- 
tion of the probe delay for three different values of the 
change in the index of refraction in the strained layer, 5n. 
As expected, the greater the change in index of refraction, 
the greater the amplitude of the differential reflectivity 
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FIG. 10: The differential reflection calculated numerically 
using Eqn. I|14[l for three different values in the change in 
index of refraction. The parameters used are = 2.65, 
uj — 3.29 eV , d = 30 nm, L — 1 /im, and C s = 7000 m/s. 
[Note that the curves are sifted to avoid overlapping.] 



oscillations. A larger change in the index of refraction im- 
plies that more electron-hole pairs are excited near the 
air-GaN interface, which acts as a stronger source for the 
coherent acoustic phonon reflectivity oscillations. These 
results are qualitatively consistent with the experimental 
results shown in Fig. In Fig. ^2 we fix the change in 
the index of refraction, Sn = 0.01, and varied the thick- 
ness of the strained layer, d. The result shows a larger 
amplitude for the differential reflection in wider strained 
layers, which is consistent with what was experimentally 
observed in Fig. |3J 

Both Fig. and Fig. suggest that the amplitude 
of the differential reflectivity seems to increase monotoni- 
cally with Sn and d. All these features can be understood 
more easily using the single-reflection approximation to 
the full formula for the reflectivity in Eqn. <39|) . As- 
suming the change in the index of refraction is small, 
\5h\ <C 1, we may select contributions from only terms 
proportional to Sh in the infinite Fabry-Perot series for 
the total reflection amplitude. The relevant reflection 
processes selected are schematically depicted in Fig. ^] 
In this case the total reflection amplitude is given by the 
leading terms in the Fabry-Perot series 



r = r Q + n + r 2 + O (Sh 2 ) 



(40) 



where ro is the background reflection amplitude, and the 
first order terms in Sh are 



„2ikz 



on, r 2 



oce l2k(z+d *>Sh. 



(41) 



FIG. 11: Calculated differential reflection for different values 
of the thickness of the strained GaN layer. The parameters 
are rib = 2.65, ui = 3.29 eV, L — 1 fim, and C s = 7000 m/s. 
The change in the index of refraction in the strained GaN 
layer was assumed to be Sn — 0.01. 





n 


ir 


n 

Cs^ 























z z+d 



FIG. 12: Schematic diagram of the single-reflection approxi- 
mation in the Fabry-Perot reflection; where we have selected 
the processes only proportional to Sh (= n — n). 



Eqn. I|14[l becomes 

= — ^ — - [Snsm(2kz + kd) + KCOs(2fcz + kd)] 

R n z — 1 

8 sin(fcd) 



n 2 — 1 



-\Sh\ sin(fcc?) sin(2fcz + cf>), 



where the phase <j) is given by 

6 = kd+ arctan ( — ) 
\5nJ 



(42) 



(43) 



Note that for small values of kd, we can expand Eq. Q42|) 
in a Taylor series to obtain 



To linear order in Sh, the differential reflectivity 



AR _ 8 kd 



n z - 1 



[Sn sm(2kz) + k cos(2fcz)] (44) 
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to first order in kd. 

From Eqn. one can see that the amplitude of the 
oscillation scales linearly with the change of the real part 
and imaginary part of the index of refraction respectively. 

Since the strain-pulse-front moves with the speed of 
sound, z = Cot, one can rewrite Eqn. Ij42(l as 

(u,t) oc \Sh\ sin(fcd) sin (^^j; + 4^J (45) 

where T = ttc/ (nbCou) — A/(2Co7ib). Note that in the 
limit kd — > 0, Eqn. (gSJ is equivalent to Eqn. ijMjl of the 
microscopic theory described in the last section. 

Equation (|45|l shows that the period of oscillations 
in the differential reflectivity is given by T and hence 
is determined by the wavelength of the probe pulse. 
This is consistent with what was seen in these experi- 
ments, but differs from previous experiments where the 
frequency of coherent phonon oscillation was determined 
by the sample and not the external probe. Note also 
that Eq. H45|l explains the experimental dependence of 
the the amplitude of the oscillation. It is proportional 
to the change in the index of refraction and varies with 
sin(/cd), and hence will increase linearly with the thick- 
ness of the strained layer as well as the probe frequency 
ui for kd = nb(ui/c)d -C 1. 

The change of index of refraction in eqn. I|45|l comes 
from two effects. The first one, which is small, is the band 
gap change (renormalization) due to the deformational 
potential phonon interaction. The other effect results 
from the change in the built-in piezoelectric field which 
leads to non-zero absorption below the band gap because 
of the Franz-Keldysh effect as shown in figure IT51 

In general, it is difficult to determine the Franz- 
Keldysh field experimentally. Here we provide a way to 
estimate the order-of-magnitude of the built-in piezoelec- 
tric field from the strength of the measured differential 
reflectivity. In Fig. 1131 we display the bulk absorption 
coefficient as a function of photon energy Tilo, in bulk 
GaN (dotted line) at a fixed piezoelectric field, F = 
0.93 MV/cm (solid line). The inset shows the absorp- 
tion coefficient at a fixed photon energy Tiui = 3.29 eV 
as a function of the piezoelectric field, F. The change 
in the absorption at Juo = 3.29 eV (below the band-gap 
Eg = 3.43 eV), can be used to estimate the piezoelectric 
field. For a probe energy below band-gap, we read off 
the amplitude of the differential reflectivity from the ex- 
perimental data. From this value, and an estimate for d 
based on information about the sample geometry, we ob- 
tain the change in the index of refraction which, in turn, 
gives an estimate for the absorption coefficient a(uj) at 
the probe photon energy In the Franz-Keldysh effect, 
the electro-absorption a(oj) is related to the piezoelectric 
field, F, viaia 

lw / [4 f (E g - huj)\ 3/2 } 

a{Lu) ^- h -E-^ e ^\-3f{^r) )■ 

(46) 



Here / = eF/(E al) where E = 3.435 eV and a = 
39.68 A are the excitonic energy and length scales in 
GaN. From this formula, for instance, we estimate that 
a piezoelectric field on the order of F « 0.93 MeV/cm 
is responsible for a 1% change in the index of refraction 
near the band-edge. 




co(eV) 

FIG. 13: The absorption coefficient as a function of the probe 
energy. The absorption tail below the band gap E g = 3.43 eV 
is due to the built-in Franz-Keldysh field. The dotted curve 
is the free-carrier absorption. The inset shows the absorption 
coefficient as a function of the piezoelectric field. The param- 
eters used in this model calculation are: d=30 nm, ^(,=2.65, 
ao~4 nm, and Eo = 15 meV. 

The corresponding changes in the real and imaginary 
parts of the dielectric function are shown in Fig. ITU 

Finally, we would like to mention the possibility of us- 
ing multiple pump pulses to coherently control the gen- 
erated coherent phonons in the epilayers. There have al- 
ready been quite a few coherent control experiments with 
multiple quantum wells from several different groups such 
as SunSfi, 6zgiis2I, and Nelson^. In the InGaN/GaN epi- 
layer systems, the effect of multiple pump pulses should 
be more prominent then in AlGaAs/GaAs systems be- 
cause of the stronger piezoelectric field. This might lead 
to novel applications in imaging or possible new devices 
based on coherent acoustic phonons. 



V. CONCLUSION 

We have presented a theory for the detection of a 
new class of large amplitude, coherent acoustic phonon 
wavepackets in femtosecond pump-probe experiments on 
In^Gai-^N/GaN epilayers and heterostructures. The 
InGaN/GaN structures are highly strained and at high 
In concentrations have large built in piezoelectric fields 
which account for the large amplitude of the observed re- 
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FIG. 14: The change in the real and imaginary parts of the 
dielectric function in GaN as a function of the probe energy 
for the situation described in Fig. 1131 The vertical dotted 
line at 3.43 eV is the band gap energy. Just below the band 
gap the change in the dielectric function is dominated by the 
imaginary part as indicated by the dotted oval. 



flectivity oscillations. This new class of coherent acous- 
tic phonons is generated near the surface and propagates 
into the structure. The frequency of the reflectivity os- 
cillations is found to be proportional to the frequency of 
the probe. These coherent phonon wavepackets can be 
used as a powerful probe of the structure of the sample. 
We are able to model the generation and propagation of 
these acoustic phonon wavepackets using a simple string 
model which is derived from a microscopic model for the 
photogeneration and propagation of coherent acoustic 
phonon wavepackets in InGaN/GaN multiple quantum 
wells. Our model successfully predicts the observed de- 
pendence of the coherent phonon reflectivity oscillations 
on probe wavelength and epilayer thickness. 
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